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In quantum dots made from materials with nonzero nuclear spins, hyperfine coupling creates a 
fluctuating effective Zeeman field (Overhauser field) felt by electrons, which can be a dominant source 
of spin qubit decoherence. We characterize the spectral properties of the fluctuating Overhauser 
field in a GaAs double quantum dot by measuring correlation functions and power spectra of the 
rate of singlet-triplet mixing of two separated electrons. Away from zero field, spectral weight is 
concentrated below 10 Hz, with ~ l// 2 dependence on frequency, /. This is consistent with a 
model of nuclear spin diffusion, and indicates that decoherence can be largely suppressed by echo 
techniques. 



Electron spins in quantum dots are an attractive can- 
didate for quantum bits (qubits) [U [2] . For gate-defined 
devices made using GaAs, the coupling of single electron 
spins to ~ 10 6 thermally excited nuclear spins creates a 
fluctuating effective Zeeman field (the Overhauser field), 
B lmc , with rms amplitude B nuc ~ 1-3 mT [3j HJ [5j [6j [7] . 
At experimentally accessible temperatures, B„„ c fluctu- 
ates both as a function of position and time, with tempo- 
ral correlations over a broad range of time scales, and is 
a dominant source of spin dephasing [8, 9- flO l fTT | [T2 j [13] 
and low-field spin relaxation [31 SI [3 HH [T5J [TD] in these 
systems. Spin manipulation schemes [IDJ H7J [TH] EFJ to 
control spin dephasing, such as spin echo and its general- 
izations, depend critically on a knowledge of correlations 
and time scales of the fluctuating nuclear environment. 

Previously, fluctuating Overhauser fields have been in- 
vestigated in atomic systems [2DJ using optical Faraday 
rotation, superconducting quantum interference devices 
[21] and force-detected magnetic resonance [22] . In quan- 
tum dots, dynamic nuclear polarization (DNP) [23J HH 
1251 126] can drive the nuclear system beyond equilibrium 
to produce fluctuating currents and feedback effects in 
connection with Pauli spin-blockade [THl [27] [25] [25] . 

In this Letter, we report measurements of the temporal 
correlations and power spectral densities of the nuclear 
environment in a two-electron GaAs double-quantum-dot 
system. In contrast to previous work JTS] [25] [23], we 
do not drive the nuclear system using DNP, but rather 
probe the statistical fluctuations of the unpolarized nu- 
clear bath in thermal equilibrium [8] [12]. Fluctuations of 
the Overhauser field are detected as fluctuations in the 
dephasing time of a two-electron spin state, making use 
of high-bandwidth proximal charge sensing |30j . Fluctu- 
ations are found to be broadband over the measurement 
bandwidth, 40 mHz to 1 kHz, and sensitive to an applied 
magnetic field in the range B = to 20 mT. Experi- 
mental results are shown to be consistent with a simple 
diffusion model of nuclear dynamics, also presented here. 

The double quantum dot is formed by Ti/Au top 
gates on a GaAs7Alo.3Gao.7As heterostructure with 



a two-dimensional electron gas (2DEG) with density 
2xl0 15 m~ 2 and mobility 20 m 2 /Vs, 100 nm below the 
surface (Fig. 1(a), inset), similar to devices reported pre- 
viously [DJ [23]. Measurements are made in a dilution re- 
frigerator with base electron temperature of ~ 120 mK. 
The conductance Gqpc of a proximal radio frequency 
quantum point contact (rf-QPC) is sensitive to the charge 
configuration of the double dot. Gqpc controls the qual- 
ity factor of an rf tank circuit, modulating the reflected 
power of a 220 MHz carrier. Demodulation yields a 
voltage Vrf, proportional to Gqpc, that constitutes the 
charge-sensing signal [3D]. The applied field, B, is ori- 
ented perpendicular to the 2DEG. 

Figure 1(a) shows the relevant energy levels of the dou- 
ble dot in the vicinity of the (2,0)-(l,l) charge transi- 
tion [first (second) index is the charge in the left (right) 
dot]. Interdot tunneling, t c , and detuning, e, from the 
charge degeneracy are determined by electrostatic gates. 
A gate-pulse (Fig. 1(b)) cycle prepares new singlets each 
iteration by configuring the device deep in (2,0), at point 
(P), where transitions to the ground state singlet, (2,0)S, 
occur rapidly |14j . Electrons are then separated to po- 
sition S in (1,1) for a time Tg where precession between 
the initial singlet and one of the triplet states is driven 
by components of the difference in Overhauser fields in 
the left and right dots, AB nuc = B l nuc - W nuc [9] H3J. 

In an applied field, the position of the separation point 
determines whether the (1,1) singlet (S) is nearly degen- 
erate with one of the (1,1) triplets, with which it can then 
rapidly mix. Mixing of S with To (the m s = triplet) 
occurs at large negative e (green line in Fig. 1(b)) where 
exchange vanishes. S — T mixing is driven by compo- 
nents of AB„ UC along the total field (applied plus Over- 
hauser fields). In contrast, mixing of S with T + (the 
m s = +1 triplet), which occurs at a less negative, field- 
dependent value of e (red line in Fig. 1(b)) where Zeeman 
splitting matches exchange, is driven by components of 
AB JlMC transverse to the total field. Measuring the degree 
of evolution out of the prepared S state following sepa- 
ration, by measuring the return probability to the (2,0) 
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FIG. 1: (Color online) (a) Schematic energy diagram of 
the two-electron system. Inset: false-color SEM image of a 
double-dot with integrated rf-QPC charge sensor similar to 
the one measured (scale bar is 500 nm). (b) Gate-pulse cycle 
that is used to prepare (P) the (2,0) singlet, separate (S) into 
(1,1), either to the S-Tq degeneracy (green dashed line) or 
the S-T+ degeneracy (red dashed line), and return to (2,0) 
for measurement (M). (c) rf-QPC readout, V t t, around the 
(l,l)-(2,0) transition during application of the cyclic gate- 
pulse sequence, showing the readout triangle indicated with 
white lines (B = mT, tq — 50 ns). A background plane 
has been subtracted, (d) Vrf as in (c), but for S at the S-T+ 
degeneracy (B = 10 mT). (e) Average value of Ps(ts) at B = 
0, tm = 2 fis. Red line is a fit to the theoretical gaussian form, 
(f) Average value of Ps(tm) showing contrast dependence, ts 
= 50 ns. Red line is a fit to the exponential form (see main 
text). 



charge configuration after a certain separation time, ef- 
fectively measures these components of the Overhauser 
field difference in the two dots. Measurement is carried 
out by moving the system to position M in (2,0) for a 
time tm = 5 [is, during which only S return to (2,0) 
with appreciable probability. The spin state — triplet or 
singlet — is thereby converted to a charge state — (1,1) or 
(2,0), respectively — which is detected by the rf-QPC. 

Figures l(c, d) show the time-averaged V r f as a func- 
tion of gate voltages Vl and Vr. Once calibrated, V r f 
gives the probability 1 — Ps that a prepared singlet 
evolved into a triplet during the separation time t$. In- 
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FIG. 2: (Color online) (a) rf-QPC sensor output Kf as a func- 
tion of Vl and Vr with gate-pulse cycle applied (ts = 25 ns, 
tm = 1.6 ps, B = 100 mT). Color scale as in Fig. 1. (b) 
Repeated slices of Vl with Vr = -709 mV as a function of 
time. Markers on left axis correspond to markers in (a), (c) 
Sensor output calibrated to Ps (blue) along with a measure- 
ment of the background QPC noise (pink) from (b) at arrow 
positions. Bandwidth limited to ~ 3 Hz. (d) Similar to (b) 
but for B = 0, color scale same as in Fig. 1. (e) Similar to 
(b) but with S-point at S-T+ degeneracy, B — 100 mT, color 
scale same as in Fig. 1. 



side the readout triangle (see Fig. 1(c)), triplet states 
remain blocked in (1,1) for a time T\ ^> tm [13] . Simi- 
larly, inside the rectangular region indicated in Fig. 1(d), 
the prepared singlet mixes with T + and becomes blocked 
in (1,1). Calibration of Vrf uses the signal in (2,0) outside 
the readout triangle, where fast, spin-independent relax- 
ation occurs via (1,0) or (2,1), to define Ps — 1, and the 
region within (1,1) to define Ps = 0. 

Fitting Ps( T s) averaged over tens of seconds 
with a gaussian [5J [T3] (Fig. 1(e)) gives T 2 * = 
h/ (g^BBnuc) ~15 ns corresponding to B nuc ~ 1.6 mT 
(N ~ 6 x 10 6 ), where g ~ —0.4 is the electron g-factor 
and /j,b is the Bohr magneton. The effect of finite T\ on 
the calibration of Pg can be accounted by introducing a 
factor C = (l-e-™/ Tl )Ti/T M [H] that relates P s to the 
value P' s corresponding to infinite T\, 1—Pg = (l — P' s )C. 
The dependence of Ps on tm (for a fixed T\ ~ 16 ys, and 
Tg = 50 ns) is shown in Fig. 1(f). Applying this factor to 
Fig. 1(e) gives P' s (ts > T 2 *) = 1/3, the expected value 
[13 , without normalizing the sensor output. 

With less averaging, P5 shows fluctuations that reflect 
fluctuations of Overhauser field components. Figure 2 
shows a slice through the readout triangle, obtained by 
rastering Vl at fixed Vr, with B — 100 mT, ts = 25 ns. 
At B = 100 mT, fluctuations in Pg have a flickering ap- 
pearance with broadband time dependence extending to 



3 



(a) 10-1 




1000 



100 



FIG. 3: (Color online) (a) Power spectra of Ps at various 
magnetic fields, rs = 25 ns. Spectra obtained by FFT (with 
Hamming window) of average of 8 traces sampled at 10 kHz. 
Background measurement noise (BG) found by setting rs = 1 
ns at B = 100 mT. Inset: numerical simulation results for cor- 
responding magnetic fields: B = (pink), B — 5 mT (blue), 
B = 10 mT (green), B = 100 mT (red), (b) Autocorrela- 
tion Ps for rs = 25 ns and B = 100 mT (red curve). Model 
function (Eq. 1) (brown) and Monte Carlo result (black). 

several seconds. Comparing the quieter (pink) trace in 
Fig. 2c, for point M such that (1,1) always returns to 
(0,2), to the fluctuating (blue) trace, where return to 
(0, 2) requires S — To mixing by Overhauser fields, we 
see that the amplitude of the fluctuating signal (blue) 
is ~ 100 times larger than the background noise of the 
charge sensor. At B = 0, slices across the readout tri- 
angle does not show a flickering (large, low- frequency) 
Ps signal (Fig. 2(d)). Figure 2(e) shows slices across the 
S — T + resonance (see Fig. 1(d)). Here also, Ps also 
does not have a flickering appearance, independent of B, 
reflecting rapid fluctuations of transverse components of 
AB nuc . We avoid rapidly cycling through the S — X+ 
transition, which can produce DNP [25] . 

To investigate the spectral content of Ps fluctuations, 
fast Fourier transforms (FFTs) of V r f are taken with Vl 
and Vr positioned to sample the center of the readout tri- 
angle. Figure 3(a) shows power spectra of Ps, with r s = 
25 ns, over the range B = - 100 mT. Measurement at rs 
= 1 ns, where Ps ~ 1, has a 1// form and is identical to 
the noise measured outside the readout triangle, and con- 
stitutes our background of instrumental noise. At B = 
no spectral content above the 1// background noise is 
seen (Fig. 2(a)). With increasing B, an increasing spec- 
tral content is observed below ~ 100 Hz. For B > 20 
mT, the spectra become independent of B. The depen- 
dence of the power spectrum of Ps on separation time rs 



is shown in Fig. 4. We found that the largest fluctuations 
over the greatest frequency range occur for rs ~ T 2 * ~ 
15 ns, and these fluctuations show a roughly l// 2 spec- 
trum. Spectra were also obtained out to 100 kHz (not 
shown) where no additional high frequency components 
were observed above the background noise. For rs < T 2 *, 
Ps remains near unity with few fluctuations; For rs > T% 
low-frequency content is suppressed while components in 
the range 1 — 10 Hz are enhanced. 

We model fluctuations in Ps as arising from the dy- 
namic Overhauser magnetic field in thermal equilibrium. 
A classical Langevin equation is used to describe fluctu- 
ations of AB„„ C arising from nuclear spin diffusion on 
distances much larger than the lattice spacing and times 
much longer than the time-scale set by nuclear dipole- 
dipole interaction. For B ^> B nuc , correlations of the 
Overhauser field can be evaluated analytically in terms of 
a dimensionless operator AP Z for each nuclear spin species 
(3, where X^^^f; = ^nuczl B nuc and similarly for 
the right dot, with x 75As = l,x 69Ga = 0.6,x 7lAs = 0.4. 
This gives (if (t + At)A£(t)) = [(1 + AtlV^) 1/2 (l + 
AtDp/a 2 ^)]^ 1 , at time difference At, where Dp is the 
species-dependent spin diffusion coefficient, a z is the elec- 
tron wave function spatial extent perpendicular to the 
2DEG (and along the external field) and a± is the wave 
function extent in the plane of the 2DEG, assumed sym- 
metric in the plane. Brackets (. . .) denote averaging over 
t and nuclear ensembles. 

Statistics of Ps for S — To mixing are found using 
the z-component of the Overhauser operators, AA Z — 
J2/3 x>3 (^z i ~ ^z,r)- For gaussian fluctuations and a 
species-independent diffusion constant, D, this gives a 
mean (Ps) — |[1 + e~ 2G ( A/1 ^] and autocorrelation 
(Ps(t + At)P s (t))-(P s ) 2 



e -iG z (AAi) 

4 



cosh(4G 2 (AA z (t- 



Ai)AA,(i))) - 1 , 

where G = ts/T^ is a gain coefficient. The autocor- 
relation function at B = 100 mT shown in Fig. 3(b) 
is obtained by Fourier transforming the power spectrum 
[3"T] . We fit to the autocorrelation function using a con- 
trast factor, C, (see Fig. 1(f) and discussion), and the 
diffusion coefficient, D, as fitting parameters. Wavefunc- 
tion widths are taken from numerical simulations of the 
device [32], a z = 7.5 nm, <r± = 40 nm. The fit gives 
D ~ 10~ 13 cm 2 /s, consistent with previous measure- 
ments on bulk GaAs samples using optical techniques 
[BJ. In Eq. (1) the dependence on rg leads to a scaling of 
the correlation time of Ps by G 2 to find the underlying 
Overhauser field correlation time. For fields B > 20 mT, 
the data in Fig. 3(b) indicate an autocorrelation time of 
~ 3 s for Ps corresponding to a time r^ ~ 10 s for AA Z 
to decorrelate by half of its initial value. 

Near B ~ 0, transverse components of the nuclear field 
lead to rapid dephasing of nuclear spins. In this regime, 
we use a Monte Carlo method to simulate nuclear dynam- 
ics |33j . Figure 3(b) shows that numerical and analytical 



4 




Frequency [Hz] 

FIG. 4: (Color online) Power spectra of P s at B = 100 mT 
for separation times rs = 25 ns (red) and rs = 100 ns (blue). 
Setting rs = 1 ns (black) yields background noise. Inset shows 
simulation results for B — 100 mT, rs = 25 ns (red) and rs = 
100 ns (blue). Note the suppression of low- frequency content 
and enhancement of mid-frequency content for long rs in the 
experiment and simulation. 

approaches agree at higher fields, where both methods 
are applicable. Numerical power spectra for B ~ are 
shown in the inset of Fig. 4. 

Experiment and theory show reduced low-frequency 
spectral content as B decreases toward zero. This can 
be understood as arising from the influence of the trans- 
verse nuclear fields at low B, which rapidly dephase nu- 
clear spins and suppress long time correlations in AB nuc . 
Similar behavior, though independent of B, is observed 
in the spectra of Ps at the S — T + resonance (not shown). 
Below B ~ 10 mT, an increased spectral content at fre- 



quencies between 1-10 Hz is observed in the experiment 
and theory. The frequency at which the spectra intersect 
however, remains constant (~ 1 Hz) in the simulations 
but increases at low B in the experimental data. We 
are able to approximate this behavior in the simulation 
by increasing the diffusion coefficient (D ~ 10~ 12 cm 2 /s 
at B = 0), implying an enhancement of diffusion, be- 
yond typical values 6 , as B approaches zero. This may 
be due to the growing influence of non-secular terms in 
the dipole-dipole interaction at low magnetic field [11 [23] . 
Diffusion maybe further enhanced at low B as a result 
of electron mediated flip-flop of nuclear spins [TSJ [M], an 
effect neglected in the simulation. 

Finally, we model how the separation time for the two- 
electron spin state affects the power spectra. Simulated 
spectra are shown in the inset of Fig. 4 for rg = 1 ns, 
25 ns and 100 ns at B — 100 mT. Good agreement with 
experiment is achieved when again accounting for the 
additional 1// noise and contrast reduction. We find that 
rs acts to filter fluctuations in AB nuc , so that for rs ^> 
T£ , low frequency correlations in AB nuc are suppressed 
in the spectra of Ps (see Eq. 1). This filtering effect leads 
to the turn-over at ~ 2 Hz evident in the spectra for r§ = 
100 ns. For ts ~ , little filtering occurs and the power 
spectra of Ps reflect the underlying intrinsic fluctuations 
of the Overhauser magnetic field. 
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